The impact of hierarchy upon the values of neutrino mixing parameters 
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A neutrino-oscillation analysis is performed of the more finely binned Super-K atmospheric, MI- 
NOS, and CHOOZ data in order to examine the impact of neutrino hierarchy in this data set 
upon the value of 613 and the deviation of 623 from maximal mixing. Exact oscillation probabil- 
ities are used, thus incorporating all powers of 613 and e := 623 — tt/4. The extracted oscillation 
parameters are found to be dependent on the hierarchy, particularly for 613. We find at 90% 



CL are A32 = 2.44+^;^;^ and 2.48+°;^^ x 10-3eV^ e = 6*23 - tt/4 = 0.06+°;?^ and 0.06+°;??, and 



-0.07 



+0.18 



and — O.I3I015, for the normal and inverted hierarchy respectively. The inverted 



hierarchy is preferred at a statistically insignificant level of 0.3 a. 
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The field of neutrino oscillations has progressed rapidly 
over the past fifteen years. The data can largely be un- 
derstood by the oscillation of the three known neutrinos 
[l|, |2| . Oscillation phenomenology invokes a unitary ma- 
trix that relates the flavor basis, in which the neutrinos 
are created or destroyed, to the mass basis, in which the 
neutrinos propagate through vacuum. This matrix can 
be written in terms of three real mixing angles, 612, O23, 
and ^13, and a phase 6 that determines CP violation. Os- 
cillations also require nonzero neutrino mass differences, 
being dependent upon the difference of the square of the 
masses, Aij := mf — m3 with rui the mass of neutrino 
i. A recent analysis Q reports the present knowledge of 
the values of the oscillation parameters. In particular, we 
note that only an upper limit exists on the size of the "re- 
actor" mixing angle |6'i3| = O.lQlg ijg (Icr), a constraint 
which arises, in part, from the null oscillation result of 
the CHOOZ experiment Q. A nonzero value of ^13 is 
requisite for the existence of CP violation in neutrino os- 
cillations; hence, we presently have no knowledge of the 
value of S. We also have no knowledge of the ordering 
of the neutrino mass eigenstate; namely, is 1113 greater or 
less than mi and m2? The former (latter) situation is 
termed the normal (inverted) hierarchy. 

Recently, hints of nonzero ^13 have been reported as a 
means by which to ease the tension between the determi- 
nation of the oscillation parameters 9i2 and A21 by the 
solar and KamLAND experiments [a, Q . At a smaller 
significance, analyses of atmospheric data also hint at 
nonzero ^13 [3, Q , though one study shows that the sig- 
nificance of the results are dependent upon the precise 
nature of the statistical analysis Q. Furthermore, an 
analysis of the updated atmospheric data finds no prefer- 
ence for nonzero ^13 [10]; however, this analysis employs 
approximate oscillation formulae and cannot be directly 
compared to these other works. Should these various 
hints be cleanly confirmed in current neutrino oscillation 
experiments [llHLSj . then one might be able to attack 
the issue of CP violation in the lepton sector. The situa- 
tion, however, is confounded by the existence of various 



degeneracies amongst the parameters [Ij, [T^ which can 
render ambiguous their extraction from experiments. To 
break the degeneracies, one needs to combine the results 
of experiments, e.g., superbeams, operating at different 
energies and /or baselines flil. Il6l - l27| . In the near future, 
the T2K [2^ and NOj^A 12911 exp eriments may be able 
to determine the hierarchy |30l - [33| . Atmospheric neutri- 
nos experiments can certainly be of some use in unravel- 
ing these unknowns as upgoing neutrinos travel through 
the earth over a large range of baselines and energies 
[25I. |34| - |42| . In Ref. [431], with fixed solar parameters, we 
analyzed the determination of the remaining oscillation 
parameters from the atmospheric, MINOS [441, CHOOZ 
[J|, and K2K [4^ experiments assuming CP is conserved 
(i.e., taking S = 0,7r); this work, unlike others, included 
the more finely binned data from Super-K [4^. With 
that data, we now analyze the effect of hierarchy on our 
parameter extraction and explore any preference for a 
particular hierarchy by the data. 

From Ref. [43, we note that vacuum oscillations are 
invariant under the map 



9i2 ^ 7r/2 - 

?13 ^^ —^13 , 
?23 ^T" 6-23 , 



'12 



A21 >-4 A21 , 
A32 M- -A31 
A31 M- -A32 



(1) 



In the limit of vanishing ^13, we see that both hierar- 
chies produce identical oscillation probabilities, provided 
the proper adjustment is made to ^12. Of course, matter 
effects [4^, |49| in the solar sector break the symmetry be- 
tween these hierarchies requiring ^12 < 7r/4, provided one 
assumes A21 > 0. To remain consistent with the solar 
data, a broken symmetry is introduced where one does 
not implement the transformation of 6*12 or ^13 in Eq. ([1]), 
but does make the transformation of the mass-squared 
differences. Given that the octant of 6*12 is known from 
matter effects, it is, in principle, possible to distinguish 
hierarchy through the precision measurement of vacuum 
oscillation channels [50 53J; however, for 6*13 ~ 0, hierar- 
chy is difficult to discern in part because of the large sep- 
aration between mass-squared differences, A21 ^ IA32I. 
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FIG. 1: (color online) The oscillation probabilities Vee and 
Ve/j. versus neutrino energy for bin I ( — 1.0 < cos 93 < —0.8) 
of the Super-K experiment using a two-density model of the 
earth. The (blue) solid curves are for both hierarchies and 
for ^13 = 0; the (red) dashed curves are for the NH and 
^13 — +0.15; the (green) dot-dot-dash curves are for the IH 
and 013 — -1-0.15. For Ve^i, the (turquoise) dot-dash curve is 
the NH and ^13 = —0.15; the (orange) dotted curve is the IH 
and 6I13 = -0.15. 



If ^13 is sufficiently nonzero, then matter effects provide 
the most promising avenue by which one might determine 
neutrino hierarchy. 

For nonzero ^13, resonant enhancement of the os- 
cillation probability "Pg/j can occur over long baselines 
which traverse the earth's mantle and/or core for ener- 
gies around 3 to 7 GeV I54i]. Using a two density model 
of the Earth [5^| given by a mantle density of 4.5 g/cm'^ 
and a core density of 11.5 g/cm"^ with radius 3486 km, 
this enhancement is apparent in Fig. [T] which shows the 
Vee and Ve^ oscillation channels for neutrinos traveling 
roughly the entire diameter of the earth. The resonance 
at 6 GeV (2.5 GeV) arises from the mantle (core) den- 
sity. Resonances only occur for neutrinos in the nor- 
mal hierarchy (NH); however, a similar resonance occurs 
for anti- neutrinos in the inverted hierarchy (IH). If 6*13 
were significantly nonzero, then one could use the pres- 
ence or absence of the resonance to discern which hier- 
archy is realized by nature. To do so, one would either 
need to have a relatively pure neutrino (or antineutrino) 
source or a detector, such as a magnetized iron calorime- 
ter, which could distinguish neutrino from antineutrino 
[llimii,!!!],!!!!. Unfortunately, the atmospheric neu- 
trino and anti-neutrino spectra are roughly equivalent, 
and water Cerenkov detectors, like Super-K, cannot re- 
solve neutrino from anti-neutrino. However, there is some 
sensitivity to hierarchy as the neutrino cross-section is 
around a factor of two greater than the anti-neutrino 
cross section over the relevant energy range. 

In Refs. [a, |43|j '^^ found that atmospheric neutrino 
data had an interesting impact on the extraction of ^13 
from the global data set. Treating the allowed parameter 



range of 613 as a continuous set (i.e., allowing negative 
values for the mixing angle is equivalent [5^, [53] to setting 
(5 = tt), it was shown that the atmospheric data placed a 
stringent upper bound on ^13. The atmospheric data only 
weakly bounds 6113 from below with the CHOOZ data 
providing the dominant constraint. As such, the data 
shows a statistically insignificant preference for 6*13 < 0. 
These constraints were traced to an excess of e-like events 
in the sub-GeV data set, a region where the oscillations 
due to the "solar" mass-squared difference is no longer 
trivial over very long baselines. Analytical expressions 
from previous studies of atmospheric neutrinos at these 
baselines and energies show that mass hierarchy will have 
at most a 10% impact upon the effective value of 6*13 in 
matter [HElHil. This is apparent in Fig.m Below 1 GeV, 
for both Pee and Pefj., the NH, (red) dashed curve, and 
IH, (green) dot-dot-dash curve, for ^13 = +0.15 nearly 
overlap. Vee is a function of 9f^ and thus the 013 == 
—0.15 curve is identical to the +0.15 curve. For "Pe^^ and 
6*13 = —0.15, the two hierarchies, the NH (turquoise) 
dot-dash curve and the IH (orange) dotted curve, also 
yield a nearly identical oscillation probability. Hierarchy 
should have little impact on the asymmetric nature of 
the bounds on ^13, as this originates from the low energy 
data. 

From 1 GeV up to 20 GeV, we see the resonance en- 
hancement of oscillations for the NH and nonzero ^13, 
the (red) dashed curve, 613 = +0.15, and the (turquoise) 
dot-dash curve, 6'i3 = —0.15. The resonance is absent for 
the IH for ^13 = +0.15, the (green) dot-dot-dash curve, 
and 6*13 = —0.15, the (orange) dotted curve. Also, we see 
that the linear in ^13 terms are small in this region. Thus 
the resonances provide information about the hierarchy 
and the magnitude of ^13, but not its sign. 

Turning to the data, in Fig. [21 we plot Ax^ versus 
the three parameters we vary: A32, e := O23 — 7r/4, and 
6*13. Here, we set Ax^ := X^ ~ Xm™ with Xm™ the min- 
imum value of x^ for the IH. The analysis tools used 
here are described in detail in Ref. [43]. This analy- 
sis incorporates the exact oscillation probabilities, and 
it makes use of the more finely-binned atmospheric data 
[3. We also include the MINOS ^ and CHOOZ [3 
data. We fix the solar parameters at their best fit values 
Q: A21 = 7.65 X 10-^ eV^ and 6*12 = 0.584. We are 
able to do this because, as can be inferred from Fig. [TJ 
the solar data is at sufficiently low energies that it makes 
no distinction between the hierarchies. To obtain the x^ 
curves we minimize the remaining varied parameters. 

Since MINOS provides the strongest constraints on 
A32 and it does not have significant matter effects, we 
expect small differences between the hierarchies, which 
we find. The minimum value of IA32I for the IH occurs 
at a slightly higher value than for the NH. The result- 
ing values at the 90% confidence level (Ax^ = 6.25 for a 
three parameter ffi) for IA32I are IA32I = 2.44^0^0 eV^ 
and 2.48+022 x 10"^ eV^ for the NH and IH, respectively. 
The fact that | A32 1 in the IH is larger than that in the NH 
is not surprising as suggested by the mappings between 
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FIG. 2: (color online) Ax^ versus the three varied param- 
eters for an analysis that utilizes the Super-K atmospheric, 
MINOS, and CHOOZ data. The solid (blue) curve is for the 
NH and the dashed (red) curve is for the IH. 



the mass-squared differences in Eq. ((T|). The difference 
in x^ between the minima for the two hierarchies is 0.9 
which is a 0.3 a effect for a three parameter fit. 

Also in Fig. [21 we present Ax^ versus e. As we vary 
ah three parameters, the difference in x^ between the 
two minima corresponding to the two hierarchies is the 
same in all three panels. The best fit values are e = 
Omto°w and OMtoTi for the NH and IH, respectively. 
This indicates a preference for 623 to lie in the second 
octant (e is positive) by Ax^ of 0.6 and 1.2, or 0.2a and 
0.4ct, for the NH and IH, respectively. 

In the last panel of Fig. [U we present Ax^ versus Oi^. 
There are differences between the hierarchies in the loca- 
tion of the minima and in the errors for ^13. Here we find 
the best-fit values for 6*13 are -O.OTtJJi^ and -O.lSta^g 
for the NH and IH, respectively. This implies that ^13 
is negative and non-zero by Ax^ of 1.8 and 2.8, or O.Scr 
and O.Sct, for the NH and IH, respectively. In Ref. |59j . 
we showed that a negative value for ^13 allows ^23 to lie 
in the second octant while still maintaining an excess of 
sub-GeV e-like events. 

These results follow what we expect from Fig. [TJ We 
find a preference for a negative 613 independent of hierar- 
chy as this arises from the excess of sub-GeV e-like events 
seen in the data; in this region, the hierarchy question is 
not relevant. In extracting this mixing angle from the 
data, the main difference between the hierarchies is that 
the minimum value of x^ shifts further from zero and the 
error bars for the IH are larger. The larger error bars 
have been previously noted \4^, J6y| . The origin of this 
and the other differences between the hierarchies arises 
from the high-energy MSW resonances present for NH 
neutrinos and IH anti- neutrinos for nonzero ^13. Due to 
the difference in cross sections, neutrinos have a greater 
impact upon the data over anti-neutrinos. As the at- 
mospheric data is less sensitive to antineutrinos, a larger 



value of [013 1 is needed for the IH to account for the data 
in the resonance region. For the normal hierarchy which 
has the resonances in the dominant neutrino channel, the 
atmospheric data is able to reasonably bound ^13; how- 
ever, for the IH, bounds from the atmospheric data on 
6*13 are less stringent. 

We have found that the present world's data have 
statistically insignificant implications for four important 
questions: the magnitude of ^13, the sign of ^13, the oc- 
tant of ^23, and neutrino hierarchy. What might be re- 
quired for atmospheric data to provide some significant 
hint of hierarchy? The answer depends on the value of 
6*13, since this controls the size of the matter resonances. 
We can provide a rough estimate by looking at the dif- 
ference between the hierarchies as predicted by the the- 
ory for each data bin for the fully contained atmospheric 
data. The upgoing neutrinos provide the greatest im- 
pact; for these data bins, we find that the theoretical 
difference between the hierarchy values is about one half 
of the present statistical error of the experimental result 
for 6*13 — ±0.15. If 6*13 is in fact this large, then a re- 
duction of the statistical error bars by a factor of two, 
or an increase in the total number of events by a factor 
of four, could begin to produce statistically significant 
indications of hierarchy. 

We have investigated the question as to how neu- 
trino hierarchy affects the extraction of ^13 from at- 
mospheric and long baseline experiments. We have in- 
cluded an exact expressions for the oscillation probabili- 
ties, which necessarily contains all linear and higher or- 
der terms in ^13 and e. We have also used the more 
finely-binned Super-K atmospheric data. The IH is pre- 
ferred at a statistically insignificant level. The extracted 
value of 6*13 differs between the two hierarchies, with 
6»i3 = -0.07;°:^? and -O.lSt";^^ for the NH and IH, 
respectively. The error on 9i3 is smaller in the NH as 
the presence of the high-energy resonances in the domi- 
nant neutrino channel leads to the Super-K atmospheric 
data restricting the value of ^13, an effect that is absent 
for neutrinos in the IH. The extracted values for ^23 and 
A32 are e = OMtoTe and OMtoT?^ ^32 = ^.U^oio and 



2.48 



-fO.25 
-0.22 



X lO'^^eV , for the NH and IH, respectively. 
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